The production of K^K pairs in proton-proton collisions below the cj) meson 

threshold 
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The pp — > ppK^K~ reaction was measured below the (p threshold at a beam energy of 2.568 GeV 
using the COSY-ANKE magnetic spectrometer. By assuming that the four-body phase space is dis- 
torted only by the product of two-body final state interactions, fits to a variety of one-dimensional 
distributions permit the evaluation of differential and total cross sections. The shapes of the dis- 
tributions in the Kp and Kpp invariant masses are reproduced only if the K~p interaction is even 
stronger than that found at higher energy. The cusp effect in the K'^ K~ distribution at the K^K'^ 
threshold is much more clear and some evidence is also found for coupling between the K~p and 
K'^u channels. However, the energy dependence of the total cross section cannot be reproduced by 
considering only a simple product of such pair-wise final state interactions. 



PACS numbers: 13.75.-n, 25.40.Ep, 13.75. Jz 



INTRODUCTION 



The original motivation for the study of kaon-pair pro- 
duction in the pp — > ppK~^K~ reaction near thresh- 
old was the investigation of the structure of the scalar 
mesons ao(980) or /o(980) Such measurements were 
initially performed by the COSY-11 collaboration at sev- 
eral different excess energies below the (/(-meson produc- 
tion threshold However, their results showed that 
scalar meson production cannot in fact be the dominant 
driving mechanism in kaon pair production and that 
the data can be explained without the explicit inclusion 
of the ao//o. Furthermore, they showed that the K^p 
and K^pp invariant mass spectra were strongly distorted, 
presumably by the K~p final state interaction (FSI) 3. 
This was most apparent in the ratio of the differential 
cross sections in terms of the K^p and A'+p invariant 
masses. 

The pp — ^ ppK^K^ reaction was also investigated 



* Electronic address: 
tElectronic address: 



qy4@phy.duke.edu] 



m . hart mann@fz-j uelich. de| 



with higher statistics above the threshold for the produc- 
tion of the (j) meson, mainly with the aim of investigating 
the properties of that meson . After removing the (p 
contribution in the spectra, it was clear that the K~p and 
K~pp distributions in the non-^ data were both strongly 
influenced by the K~p interaction d, It has been 
suggested that this is connected with the production of 
the A(1405) excited hyperon Q, which might be treated 
as a KN quasi-bound state with a width that overlaps 
the KN threshold This idea was put on a quanti- 
tative footing by assuming that the A(1405) was formed 
through the decay N* iv:+A(1405) JT,]. The strength 
and details of the KN interaction are clearly important 
elements in the interpretation of possible kaon nuclear 
systems, such as the deeply bound K'pp states [l^ . 

In addition to the K^p FSI, and one between the two 
protons, the data also showed an enhancement at low 
K^K~ invariant masses with some possible structure at 
the K^K'^ threshold ^-JJ. Though the effects are small, 
they might be influenced by the ao(980) or /o(980) scalar 
mesons. However, the investigation of this region was 
hampered by the need to separate the non-0 from the (j) 
contribution and the fact that the data were spread over 
a very wide range of K^ K~ invariant masses. Measure- 
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ments below the 4> threshold can provide useful informa- 
tion on these interesting FSI effects without suffering the 
distortion of the (j) meson. However, the limited statistics 
in the low-energy COSY-11 data are insufficient for 
detailed studies. 

Previous measurements of the pp — >■ ppK^ K~ reac- 
tion were carried out at the COSY-ANKE magnetic spec- 
trometer at e = 51, 67, and 108 MeV [I,!], where the 
threshold is at e — 32.1 MeV. Here the excess energy is 
defined as e = y/s — 2{mp + ■mK)c^, where a/s is the total 
center-of-mass energy and rup and rrifc are the particle 
masses in the final state. Because of the limited accep- 
tance of this spectrometer, an ansatz has to be made 
regarding the distribution of events over the four-body 
phase space in order to convert count rates into cross sec- 
tions. This was done assuming that the distortions were 
the products of those present in the two-particle subsys- 
tems. All the ANKE non-0 data seemed to be consistent 
with an effective scattering length of a-K-p = (0-1-1. 5i) fm 
with no obvious influence of an energy dependence associ- 
ated with an effective range term. The dominance of the 
imaginary part is not unexpected because of the strong 
couplings to the Ett and Att channels but, due to the 
presence of two other final-state particles, this parame- 
ter is not necessarily an intrinsic feature of the isolated 
K~p system. 

The ANKE measurements at three excess energies 
also showed some enhancement at low K'^K^ invari- 
ant masses but with at least a break of slope at the 
K^K^ threshold. A combined analysis of all the results 
in this region 0, [3 shows that the data can be under- 
stood in terms of a final state interaction involving both 
K^K~ elastic scattering plus a contribution from the 
K'^K~ ^ K^K^ charge exchange. Although suggestive, 
the data are not sufficient to draw firm conclusions. 

In this paper we present much more precise pp — > 
ppK~^K~ differential cross section data at a beam en- 
ergy of Tp = 2.568 GeV (e = 23.9 MeV) obtained using 
the COSY-ANKE spectrometer. With high statistics on 
the reaction below the (jj-meson threshold, we can study 
the effects of the final state interactions in the K^p and 
K'^K" systems in greater details. 

The paper is organized as follows. We first describe the 
experimental setup and data analysis in Sec [Hi Given 
that the procedures involved are similar to those em- 
ployed at higher energies d, Q , this can be quite brief. 
The fitting of the phenomenological parametrization to 
the raw pp — >■ ppK'^K~ data in order to make accep- 
tance corrections is also described here. The resulting 
differential cross sections and total cross section for the 
pp — >■ ppK~^K~ reaction are presented in Sec lIIIl followed 
by our conclusions in Sec IIVI 



II. EXPERIMENT AND DATA ANALYSIS 

The measurement of the pp — > ppK~^K~ reaction was 
performed at an internal target station of the Cooler Syn- 



chrotron (COSY) of the Forschungszentrum Jiilich [13| . 
The ANKE spectrometer [3 which consists of 

three dipole magnets, registers positively and negatively 
charged ejectiles in the side detection systems, with the 
fast positively charged particles being detected in the for- 
ward detector. Particle identification relies on time-of- 
flight measurements [1, HMtI from START and STOP 
counters, and momentum information obtained from the 
multiwire proportional chambers. 




0.92 0.94 0.96 

Missing mass (K*K"p) [GeV/c^] 

FIG. 1: The pK'^K~ missing-mass distribution in the pp — ^ 
pK+K'X reaction at Tp = 2.568 GeV. The hatched his- 
togram shows the cuts imposed for the selection of the non- 
detected proton. The solid line, which is a second-order poly- 
nomial fit, was used to estimate the background contribution 
under the proton peak. 

In close-to-threshold production experiments, the to- 
tal cross section changes very rapidly with small changes 
in the excess energy. The proton beam energy, Tp — 
2.568 GeV, was therefore determined very prec isely 
through a careful monitoring of the Schottky spectra [131 ■ 
The resulting value of the excess energy with respect to 
the ppK^K" production threshold, e — 23.9 MeV, is 
well below the nominal (j) threshold. However it should 
be noted that this is an average value, since the beam 
energy decreases by up to 4.6 MeV through the course of 
a machine cycle due to the interaction with the target. 
This effect was also investigated in the simulation. 

The experiment relied on a triple-coincidence, involv- 
ing the observation of a K^K~ pair in the side detec- 
tors and a fast proton in the forward detector. The 
pp — >■ ppK~^K~ reaction was then identified by requir- 
ing that the missing mass of the K^K^p system be con- 
sistent with that of a proton. In the analysis, a ±3tT 
((7 = 2 MeV/c^) cut was applied on the missing-mass 
distribution of the selected K^K~p events, as shown in 
Fig. [TJ The fraction of misidentified events inside the 
cut window around the proton mass was estimated to 
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be about 5%, which was subtracted from the peak using 
weighted data from the side bands, as parameterized by 
the sohd hne. Any ambiguity in this procedure, which is 
less than 3%, is one source of systematic uncertainty. 
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FIG. 2: DifFerential distributions of experimental (points) and 
simulated (histograms) yields for kaon pair production in the 
pp — >■ ppK'^ K~ reaction at e = 23.9 MeV. Vertical bars rep- 
resent the statistical uncertainties and horizontal ones the bin 
widths. The individual panels are (a) the cosine of the polar 
angle of the in the K~ reference frame, (b) the polar 
angle of the kaon pairs in the overall cm. frame, (c) the polar 
angle of the emitted proton in the pp reference frame relative 
to the beam direction, (d) the polar angle of the proton in the 
pp reference frame relative to the direction of the kaon pair, 
(e) the proton momentum in the pp reference frame, and (f) 
the K^K~p invariant mass distribution. 



After identifying clean pp — > ppK^K^ events in 
ANKE, acceptance corrections must be performed in or- 
der to evaluate differential cross sections. The simple 
ansatz used on data taken above the meson produc- 
tion threshold tried to take into account the influence of 
final state interactions in the various two-particle sub- 
systems [1, This ansatz, which is also the basis for 
the current simulation, assumes that the overall enhance- 
ment factor F is the product of enhancements in the pp, 



K~^K , and K p systems: 

F = Fpp{qpp) X FKp(qKpi) x FxpiqKpi) x FxRi^KK), 

(1) 

where qpp, qxpi, qKp2i and qKK sue the magnitudes 
of the relative momenta in the pp, the two K~p, and 
the K~^K~ system, respectively. It is believed that the 
K~^p interaction might be weakly repulsive and, if so, its 
neglect would be interpreted as extra attraction in the 
K~p system. The FSI enhancement in the K~p case 
was calculated in the scattering length approximation, 
Fxpiq) ~ 1/]! ~ iqa]"^ and the best fit to the higher 
energy data [y, Q was found with a purely imaginary 
effective scattering length, a^-p « fm. The proton- 
proton enhancement factor was derived from the Jost 
function [1, Q. The enhancement factor in the K^K~ 
system takes into account elastic K'^K~ scattering plus 
the charge-exchange K'^K~ ^ K^K^ 

The seven degrees of freedom required to describe the 
unpolarized ppK~^K~ final state were chosen to be four 
angles, the K'^K~ and K~^K~p invariant masses, and 
the relative momentum of the protons in the pp sys- 
tem Q . Distributions in these seven variables were 
generated inside the ANKE acceptance and compared 
with the experimental data, some of which are shown 
in Fig. [21 The best fit to the data was achieved with 
o-K-p — (2.45 ± OA)i fm, which is significantly larger 
in magnitude than the starting value of ax-p = l-5i fm. 
The uncertainty in the real part is large and strongly cor- 
related with the imaginary part. To allow easy compari- 
son with the analysis of the higher energy data [1, [1] , the 
effective scattering length was taken to be purely imagi- 
nary. 

The KK scattering lengths for isospin-one and zero 
were taken as in our previous work [7,] and the ratio of 
the 1=1 and / = production amplitudes of s-wave 
KK pairs was parameterized as Ce"^'=. The best fit was 
obtained with C = 0.54±0.03 and = -112°±4°, which 
are consistent with our earlier evaluation 0, [1] based on 
the above (p threshold data. The resulting descriptions 
of the experimental data in Fig. [2] are very good and cer- 
tainly sufficient for evaluating the acceptance corrections. 

The luminosity needed in the analysis was determined 
with an overall systematic uncertainty of 9% by measur- 
ing pp elastic scattering in the forward detector . This 
was checked by simultaneous studies of the beam current 
and Schottky spectra [3, which could fix the absolute 
luminosity with a systematic uncertainty of 6%. Within 
these uncertainties the two methods agreed but, in order 
to be coherent with our previous work, the luminosity 
extracted from the pp elastic scattering data was used in 
the final analysis. 



III. RESULTS 

The differential cross section for the pp — > ppK'^K~ 
reaction at an excess energy e = 23.9 MeV is shown in 
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FIG. 3: (Color online) The K invariant mass distribution 
for the pp ppK+K~ reaction at e = 23.9 MeV. The dotted 
curve shows the four-body phase space simulation whereas 
the inclusion of the final state interactions through Eq. ^ 
gives the dashed curve for flj^-^ = 1.5i fm and the red solid 
curve dK-p = 2.45i fm. The dot-dashed curve was obtained 
by considering only the pp and K~p final state interactions 
with a^-p = 2.45i fm. 
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FIG. 4: Ratio of the measured invariant mass distri- 

bution in the pp ppK^ K~ reaction to a simulation that in- 
cludes only K~p and pp final state interactions (shown by the 
dot-dashed curve in Fig. In addition to the current data 
(solid circles), weighted averages of previous measurements 
(open squares and circles) are also presented. The solid curve 
represents the best fit in a model that includes elastic K'^ K~ 
FSI and K'^K° ^ K+ R- charge-exchange 0. The best fits 
neglecting charge exchange and including only this effect are 
shown by the dashed and dot-dashed curves, respectively. 



Fig. |3] as a function of the K^K^ invariant mass. Also 
shown are simulations based on a four-body phase space 
and this distorted by the final state interactions in the 
K^K^ , pp, and K^p systems within the product ansatz 
of Eq. ([ij . This was done separately with effective scat- 
tering lengths of a^-p = l-5i fm and a^-p = 2.45i fm. 

The most striking features in the data are the strength 
near the K^K^ threshold and the dip at M^+x- ~ 
0.995 GeV/c^, which corresponds precisely to the K^K^ 
production threshold 7] . This is compelling evidence for 
a cusp effect coming from the K^K^ ^ K^K^ transi- 
tions. To investigate this phenomenon in greater detail, 
the K'^K^ invariant mass distribution was divided by a 
simulation where only the final state interactions in the 
pp and K^p, with a^-p = 2A5i fm, were considered. 
The best fit to the data shown in Fig. 0] is achieved with 
a contribution from the isospin-zero channel that is about 
three times stronger than the isospin-one. This finding is 
consistent with our earlier result Q. The deviations ap- 
parent in Figs. [3] and m at high K'^K" invariant masses 
might be connected with the approximations made in our 
coupled-channel model [7j . 

Previous analyses of the pp — )■ ppK^K~ reaction at 
different excess energies [1, [a, H have all shown a 
strong preference for low values of the K~p and K~pp 
invariant masses, Mj^~p and M^-pp- To study this fur- 
ther, we have evaluated differential cross sections in terms 



of these invariant mass distribution and also their ratios: 

_ da/dMii-p 
" da/dMK+p ' 
da/dMK-pp 

^Kpp — , (,^J 

da/dMK+pp 

The corresponding experimental data and simulations 
are shown in Figs. [5]and|ni Both Rkp and Rkpp display 
the very strong preferences for lower invariant masses 
seen in the earlier data. The low mass enhancements 
in Figs. [SJ: and IB}; clearly indicate once again that the 
pp — > ppK^ K~ reaction cannot be dominated by the 
undistorted production of a single scalar resonance oq or 
/o. Within a four-body phase space simulation both ra- 
tios should be constant and equal to one and such a simu- 
lation also fails to describe the individual Mkp and Mkpp 
distributions. Whereas the inclusion of a K~p FSI with 
an effective scattering length a^-p = 1.5i fm improves 
the situation, it overestimates the data in the high in- 
variant mass regions for both Rxp and Rkpp- With the 
larger effective scattering length a^-p — 2Abi fm, these 
ratios, as well as the individual distributions, can be well 
reproduced. Within the product ansatz of Eq. ^ the 
K~p final state interaction effectively becomes stronger 
at lower excess energies. This illustrates the limitations 
of this simple ansatz to the complex four-body dynamics. 

Although the K~p elastic final state interaction de- 
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FIG. 5: (Color online) Differential cross sections for the 
pp — > ppK^K~ reaction with respect to the invariant masses 
of K^p (upper panel) and K~p (middle panel), and their ra- 
tio Rkp (lower panel). The dotted histograms represent the 
four-body phase-space simulations, whereas the red solid and 
dashed black ones represent estimations based on Eq. lU that 
take into account K~p, pp and K~ final state interactions 
with aji-„ = 2.45i fm and a^- = 1.5i fm, respectively. 



scribes well the vast bulk of the data shown in Figs. [5] 
and [6l it is interesting to note that there seems to be 
a small but significant deviation between the K~p data 
and simulation in Fig. [SJj at low invariant masses. Since 
the threshold is at 1437.2 GeV/c^, this suggests 

that the data in this region might also be influenced by 
K~p ^ K'^n channel coupling. 

Due to the low statistics, the COSY-11 data at 10 and 
28 MeV [1, [l^ cannot distinguish between predictions 
based on effective scattering lengths of a^-„ = 1.5z fm 
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FIG. 6: (Color online) Differential cross sections for the 
pp — >■ ppK^K^ reaction with respect to the invariant masses 
of K^pp (upper panel) and K^pp (middle panel), and their 
ratio Rkpp (lower panel). The conventions for the theoretical 
estimates are as in Fig. [5] 



and ax-p = 2.4:5i fm. This illustrated for the Rkp ratio 
in Fig. [7] but this lack of sensitivity is equally true for 
Rkpp- 

The acceptance-corrected invariant mass distribution 
shown in Fig. [3] was used to determine the value of the 
total cross section, a = 6.66 ± 0.08 ± 0.67 nb, where 
the errors reflect the statistical and systematical uncer- 
tainties, respectively. The systematic effects considered 
here arise from the background subtraction, acceptance 
correction, tracking efficiency correction, and luminosity 
determination. 

The total cross section result is plotted in Fig. [S] along 
with previous measurements from DISTO Q, COSY- 
11 [Ji, [H, and ANKE d, i]. The new point seems 



6 



■o 



T3 



10' r 




1.434 1.436 1.438 

Kp invariant mass [GeV/c^] 



■D 



■o 



10' r 




1.44 1.45 

Kp invariant mass [GeV/c^] 



FIG. 7: (Color online) The distributions of the ratio Rkp for 
the pp -> ppK^K~ reaction measured by COSY-11 at (a) 
£ = 10 MeV and (b) 28 MeV [l9|, respectively. The dotted 
istograms represent the four-body phase-space simulations, 
whereas the red solid and dashed ones represent the theo- 
retical calculations taking into account K~p, pp and K'^ K~ 
final state interactions with a^-p = 2.45i fm and 1.5i fm, 
respectively. 
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FIG. 8: Total cross section for the pp — >■ ppK K~ reac- 
tion as a function of excess energy e. The present result 
(cross) is shown together with earlier experimental data taken 
from DISTO (triangle), ANKE (circles), COSY-11 (squares). 
The dotted line shows the four-body phase space simula- 
tion, whereas the solid line represent the simulations with 
ttj^-p = 1.5i fm. The predictions of a one-boson exchange 
model are represented by the dashed line [20] . 



cantly underestimated. An increase in the value of aj^-p 
might help in this region but the coincidence of strong ef- 
fects in different two- or even three-body channels must 
also bring the factorization assumption of Eq. (HI) into 
question. The dashed line, vifhich represents a calcula- 
tion within a one-boson exchange model HO], also un- 
derestimates the near-threshold data. This model in- 
cludes energy dependent input derived from fits to the 
K^p — >■ K^p total cross sections, though it does not 
include the pp final state interaction. 



high compared to the COSY-11 result at e = 28 MeV 
but it should be noted that this value was already in- 
creased by 50% compared to that originally quoted. This 
was achieved through a re-analysis of the data that in- 
cluded pp and K~p final state interactions with aj^-p ^ 
l.bi fm [19]. It is therefore possible that a further increase 
might follow if the K^K~ final state interaction were in- 
cluded in the evaluation of the acceptance corrections as 
well as an effective scattering length of a^-p — 2.45i fm 
was used instead. There seems to be less sensitivity to 
the FSI assumptions at e = 10 MeV, where the original 
re-analysis only increased the total cross section by 20%. 

It is clear from Fig. [5] that the four-body phase space 
cannot reproduce the energy dependence of the total 
cross section. With the inclusion of the pp, K^K~, 
and K~p FSI, with an effective scattering length of 
a^-p = l-5i fm, the data above the (j) threshold can 
be described well but those at lower energy are signifi- 



IV. DISCUSSION AND CONCLUSIONS 

The production of K~ pairs has been measured 
in the pp — > ppK^K^ reaction channel at an excess 
energy of e = 23.9 MeV. Even taking its 4.3 MeV/c^ 
width into account, this is well below the central 0- meson 
threshold at 32.1 MeV. The reaction was identified in 
ANKE through a triple coincidence of a K^K~ pair and 
a forward-going proton, with an additional cut on the 
K^K^p missing-mass spectrum. The high statistics and 
low excess energy allow us to produce a detailed K'^K~ 
invariant mass distribution below the cf) threshold. 

The distortion of both the K~p and K~pp spec- 
tra, which are even stronger than in our higher energy 
data, can be explained quantitatively within the product 
ansatz of Eq. ([T|) with an effective K~p scattering length 
a^-p ~ 2.45i fm. This is to be compared with the l.bi fm 
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obtained from the analysis of data measured above the (f) 
production threshold. A full treatment of the dynamics 
of the four-body ppK^K^ channel is currently imprac- 
tical. As a consequence, an energy dependence of a^-p 
is possible because this is merely an effective parameter 
within a very simplistic description of the four-body final 
state interaction. The strong K^p final state interaction 
may be connected with the A(1405) in the production 
process and it has been suggested [ll| that the produc- 
tion of non-cj) kaon pairs proceeds mainly through the as- 
sociated production pp — >■ K^pA{li05). This would also 
lead to deviations from the simple product ansatz for the 
final state interactions, not least because an attraction 
between the A(1405) and the proton would involve three 
final particles. 

Our results show a very strong preference for low K^pp 
masses and this effect seems to be even more marked than 
in the higher energy data [gI, Although this might 
be connected with the ideas of a K^pp deeply bound 
state [13, [2T| - [23j , it must be stressed that our data were 
measured far above threshold. They should not therefore 
be taken as necessarily implying that the K~ will bind 
with two protons. 

There is strong evidence for a cusp effect arising from 
the K^K'^ ^ K^K~ transitions. Our analysis within a 
coupled-channel description suggests that, with the val- 
ues of the KK scattering lengths used, the production 
of isospin-zero KK pairs dominates. Though this is con- 



sistent with results extracted from data taken above the 
(f) threshold 0, [1] , there is clearly room for some refine- 
ment in the model. On the other hand, the structure 
of the K~p invariant mass distribution of Fig. [S] in the 
1437 MeV/c^ region suggests that there might be im- 
portant coupling also between the K~p and K'^n sys- 
tems. It is evident that the interactions in the four-body 
ppK~^K~ final state are extremely complex. Neverthe- 
less, the energy dependence of the total cross section can 
be well described above the (f> threshold by introducing 
the effects of the pp, K~^K~ and K~p final state interac- 
tion with an effective scattering length of a^-p — l-5i fm. 
This would, however, have to be increased to have any 
hope of fitting the lower energy data. Further theoretical 
work is required to clarify the reaction mechanisms. 
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